Particle- based modeling of oxygen discharges 
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We present an one-dimensional particle-in-cell Monte-Carlo model for capacitively coupled 
radio-frequency discharges in oxygen. The model quantitatively describes the central part 
of the discharge. For a given voltage and pressure, it self-consistently determines the electric 
potential and the distribution functions for electrons, negatively charged atomic oxygen, and 
positively charged molecular oxygen. Previously used collision cross sections are critically 
assessed and in some cases modified. Provided associative detachment due to metastable 
oxygen molecules is included in the model, the electro-negativities in the center of the 
discharge are in excellent agreement with experiments. Due to lack of empirical data for the 
cross section of this process, we propose a simple model and discuss its limitations. 



1. Introduction 

Discharges in reactive gases such as oxygen 
play an important role in plasma-assisted etch- 
ing and thin-film deposition techniques. The re- 
quirements on the controllability and reliability of 
these discharges are so high, that further advance- 
ment of this technology critically depends on im- 
proved descriptions of the physical processes. In 
particular, the interplay between the macroscopic 
electrodynamics, which is used to control the dis- 
charge, and the microscopic plasma-chemical pro- 
cesses, which sustain the discharge and give rise 
to the materials processing, has to be understood 
not only qualitatively but quantitatively. 

Particle-based modeling is well suited for this 
task, because it directly simulates the Boltzmann- 
Poisson system describing the discharge with- 
out any assumptions concerning the species' dis- 
tribution functions or the electric field. Pro- 
vided the elementary plasma-chemical processes 
are well characterized in terms of cross sections, 
they can be easily incorporated in the collision in- 
tegral of the Boltzmann equation. Through the 
source term of the Poisson equation, the plasma- 
chemistry is then linked to the electrodynamics of 
the discharge. 

In the following we give a brief account of our 
implementation of the particle-in-cell Monte- 
Carlo collision (PIC-MCC) approach for the 
modeling of capacitively coupled radio-frequency 
(rf) discharges in oxygen. Our main focus 
will be the critical assessment of cross section 
data for (0^,02) charge exchange scattering. 



ion-ion neutralization, and detachment due to 
metastables. A full description of our approach, 
whose treatment of collisions is closely related to 
the direct simulation Monte-Carlo approach for 
rarefied gases, will be given elsewhere [1]. 

2. Model 

The description of an oxygen discharge could 
be based on a brute force numerical solution of 
the Boltzmann-Poisson system which couples the 
distribution functions of the relevant species with 
the electric potential. In most cases, however, 
this approach is not practical. More promising 
are methods which track the spatio-temporal evo- 
lution of a sample of pseudo-particles subject to 
elastic, inelastic, and reactive collisions. These 
approaches are based on the decoupling of col- 
lisions from the free flights in the self-consistent 
electric field. When the cross sections for the col- 
lisions are known, the durations of the free flights, 
as well as the probability for a collision of a par- 
ticular type to occur, can be simply obtained from 
elementary kinetic considerations. 

Our simulations of capacitively coupled rf dis- 
charges in oxygen [2,3] are restricted to the central 
axial part of the reactor (see Fig. [1]). Ignoring 
the (electric) asymmetry between the grounded 
and powered electrode, we use an one-dimensional 
(ID) model, which keeps only one spatial coordi- 
nate, < X < L, where L is the distance between 
the electrodes, but retains all three velocity co- 
ordinates Vx,Vy, and Vz- One of the electrodes 
is electrically driven by a time-dependent volt- 
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Fig. 1: Schematic geometry of the rf discharges used 
in Refs. [2,3]. We use an ID model to simulate the 
central axial part of the discharge (thick solid line). 

age, Ur{{t) = U sin27r/, while the other is set to 
U = 0. Both electrodes are totally absorbing, 
secondary electron emission is neglected, and the 
oxygen molecules are treated as an inexhaustible 
reservoir, the density of which is uq^ = p/kT, 
where T = 300 K and p is the gas pressure. 

The complex plasma-chemistry of oxygen gives 
rise to a large variety of collisions. In Table [1] 
we show the ones with the largest cross sec- 
tions. They are included in our model. We 
simulate only three species: electrons (e), nega- 
tively charged atomic oxygen (0~), and positively 
charged molecular oxygen (O^). Neutral parti- 
cles are not treated kinetically. They are only ac- 
counted for in as far as their production leads to 
an energy loss for electrons (collisions (10)-(15)) 
and in as far as they affect the balance of simu- 
lated particles (reactions (16)-(22)). 

Our collection of cross sections is semi- 
empirical, combining measured data with models 
for the low-energy asymptotic. A complete dis- 
cussion of the molecular physics entering the sim- 
ulation will be given elsewhere [1]. Here it suffices 
to mention that the cross sections for (02,0^) 
charge exchange scattering (9), ion-ion neutral- 
ization (17), and detachment (19,20) significantly 
deviate from the ones used previously [4] (see 
Fig- EI)- Our simulations indicate that the modi- 
fications are essential for obtaining results in ac- 
cordance with experiments [2,3]. 

Using acx[E) = cJm(-E')/2, where am and 
(Tex denote, respectively, the momentum and 
charge exchange cross section, and E is the rel- 
ative kinetic energy, we based our cross section 
for (©2,0^) charge exchange scattering below 
0.251 eV and above 8.5 eV on empirical data for 
momentum scattering [5,6]. For energies in be- 
tween, we employed a linear interpolation. With 
this cross section, we obtained velocity distri- 



Table 1: Elastic, inelastic, and reactive collisions in- 
cluded in our model [1]. 
elastic scattering 

(1) e + e ^ e + e 

(2) O- + 0- ^ O- + 0- 

(3) 0+ + 0+^0+ + 0+ 

(4) e + 0+ ^ e + 0+ 

(5) e + O" ^ e + 0~ 

(6) O- + 0+ ^ O- + 0+ 

(7) e + 02 ^e + 02 

(8) O- + O2 ^ O- + O2 

(9) 0+ + O2 ^ O2 + 0+ 
electron energy loss scattering 

(10) e + 02^e + 02(1^ = 1,. ..,4) 

(11) e + 02^e + 02(Ryd) 

(12) e + O2 ^ e + 0(3P) + 0(3P) (6.4 eV) 

(13) e + 02^c + 0(3P) + 0(lD) (8.6 eV) 

(14) e + 02^e + 02(aiAg) 

(15) e + 02 ^e + 02(biSg) 
electron &: ion production & loss 

(16) e + 0^ ^0 + 

(17) O- + 0+^0 + 02 

(18) e + 02 ^ + 0- 

(19) O- + O2 ^ + 02 + e 

(20) O- + 02(a^Ag) ^ 03 + e 

(21) e + 02 ^ 2e + 0^ 

(22) e + O" ^ + 2e 



bution functions in close agreement with experi- 
ments [2] . With the charge exchange cross section 
given in Ref. [4], on the other hand, we could 
not reproduce the experimental findings - nei- 
ther with our PIC-MCC code nor with the BITl 
code [7], which we used to cross-check our results. 

For ion-ion neutralization, we employed the 
cross section of a two-channel Landau-Zener 
model [8], cr„(S) = 47ri?^.(l + -^), where 
is a free parameter which we adjusted to obtain 
the experimentally measured cross section at high 
energies [9]. This cross section deviates dramati- 
cally from the one used in Ref. [4] , but it is based 
on a clear physical picture for the neutralization 
process as well as empirical data at high energies. 

Detachment of O^ on neutrals occurs through 
direct detachment (19) and associative detach- 
ment (20). The latter is rather surprising be- 
cause there is no evidence for it in beam experi- 
ments [10]. Yet, investigations of O2 discharges 
strongly suggest that associative detachment is 
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Fig. 2: Cross sections for (0^,02) charge exchange 
scattering (9), ion- ion neutraUzation (17), direct de- 
tachment (19), and associative (20) detachment. The 
dotted hnes indicate the respective cross sections used 
in Ref. [4]. 



Fig. 3: Cycle-averaged electron and ion density pro- 
files for a rf discharge in O2 without associative de- 
tachment (20) taken into account. The parameters 
of the discharge are L = 2.5 cm, / = 13.6 MHz, 
p = 13.8 Pa and U = 250 V. 



possible because of the presence of metastable 
©2(0^ Ag) and may be even the main loss process 
for 0~ in some pressure range [3]. Due to lack 
of empirical data for this process, we employed a 
simple model, which describes detachment as the 
"inverse" of a Langevin-type electron capture into 
an attractive auto-detaching state of The 
cross section of which reads 



a^aiE) = 5.96 • 
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where we assumed the polarizability of 62(0^ Ag) 
to be the same as for 02- Note, in contrast 
to the cross section for direct detachment, the 
cross section for associative detachment has no 
threshold. To determine the probability for this 
process, we also need the density of 02(a^Ag). 
Within the three species model this density is 
unknown. It should be however of the order of 
the O2 density. As a first step, we write therefore 
n/\ = C ■ no2) with C < 1 an adjustable fit 
parameter. 

3. Results 

Similar to other electro-negative gas discharges, 
the presence of negative ions in an oxygen dis- 
charge leads to abrupt changes in the ion den- 
sity which, in most cases, forces the discharge to 
stratify into a quasi-neutral ion-ion and a periph- 
eral electro-positive edge plasma. The details of 
the stratification depend on the interplay between 
plasma-chemistry and electrodynamics. 

In oxygen discharges two plasma-chemical pro- 
cesses are of particular importance: Ion-ion 



neutralization and associative detachment due 
02(a^Ag). Although the three species plasma 
model cannot fully describe associative detach- 
ment, because it assumes an homogeneous back- 
ground of 02(a^Ag) molecules, it nevertheless 
gives clear evidence that the latter process is in- 
dispensable for a correct description of experi- 
ments. 

To demonstrate this, we simulated the dis- 
charge of Katsch and coworkers [3]. In Figs. [3] 
and m we show, respectively, the quasi-stationary, 
cycle- averaged density profiles of the simulated 
charged particles for p = 13.8 Pa, U = 250 V, 
L = 2.5 cm, and / = 13.6 MHz without and 
with associative detachment taken into account. 
The parameter C ~ 1/6, implying that one-sixth 
of the O2 molecules are in the metastable state. 
The precise value of C should not be taken too se- 
rious because it is based on a rather crude model 
for associative detachment. More important is 
that without this process {C = 0), the simulation 
could not reproduce the measured densities. 

Whereas the simulation with associative de- 
tachment reproduces reasonably well the densi- 
ties of charged particles in the center of the dis- 
charge, the shapes of the (axial) density profiles 
deviate from the measured ones. Compared to ex- 
periment, the central plasma is too narrow, most 
notably, for lower voltages (not shown here) [1]. 
This is a shortcoming of the three species plasma 
model which ignores the spatial dependence of the 
02(a^Ag) density which, in reality, results from 
the interplay of volume and surface loss and gen- 
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Fig. 4: Same plot as Fig. [3] but now with associative 
detachment (20) taken into account. Experimentally, 
He ~ Uq- ^ 0.9 • 10^ cm~^ (grey bullet) [3]. 

eration processes. Because the probability of as- 
sociative detachment is proportional to riA, the 
02(a^Ag) density profile should strongly affect 
the density profiles of charged particles. 

In Fig. [5] we finally plot the time-resolved 
current densities through the discharges shown 
in Figs. [3] and m The current densities are almost 
the same, irrespective of whether associative 
detachment is taken into account or not. This 
is a consequence of the fact that in both cases 
the main part of the current is carried by the 
electrons whose densities, in turn, are basically 
identical. Further studies are needed to reveal if 
this is an artifact of our model, or if indeed two 
rather different density profiles are consistent 
with a given external power supply, (C/rf ■jrf)rf 
where A is the area of the electrodes and (...)rf 
denotes the cycle average. The configuration 
realized in the discharge depends then on the 
plasma-chemistry, in particular, on the outcome 
of the competition between ion-ion neutralization 
and associative detachment. 

4. Conclusions 

We constructed an ID three-species PIC-MCC 
model for capacitively coupled rf discharges in 
oxygen. The model simulates three kinds of 
charged particles (e, 0~, and O^) and retains 
neutral particles indirectly via collisions with the 
simulated charged particles. This is sufficient to 
reproduce measured central electron and ion den- 
sities. However, the (axial) ion density profiles of 
the simulations are too narrow compared to the 
experimental ones. We expect better agreement, 
when the modeling also allows for the possibility 
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Fig. 5: Time dependence (for five rf cycles with dura- 
tion Tj-f ) of the current densities corresponding to the 
cycle-averaged, quasi-stationary density profiles shown 
in Figs. [3] and m respectively. 

of 02(a^Ag) density profiles. For that purpose, 
the 02(o^Ag) molecules, together with their 
main loss and generation processes, have to be 
explicitly included in the Monte-Carlo collision 
approach. 
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